% Plot the power at M2 (12.421 hours) of the observatory data on to the
% predicted power map

S = dir('/home/mnair/projects/tides/*_tidal_spectra.mat'); % Only PMTM results (variable P) is useful

periods = [4 4.8 6 8 11.967236 12  12.421 12.6583 23.934472 24 25.891];
for nstn = 1:length (S);
    
    eval(['load /home/mnair/projects/tides/' S(nstn).name ]); %pmtm results
       
    station = S(nstn).name(1:3);
    
    eval(['save /home/mnair/projects/tides/' S(nstn).name(1:3) '_tidal_spectra_new latitude longitude st en P F spectra_data_ls_m2_array spectra_all_ls spectra_night_ls spectra_data_night_ls spectra_data_all_ls spectra_data_night_rob spectra_night_rob spectra_all_rob spectra_data_all_rob periods station' ]); %pmtm results
    
end;

